Dark state cooling of atoms by superfluid immersion 
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We propose and analyse a scheme to cool atoms in an optical lattice to ultra-low temperatures within a Bloch 
band, and away from commensurate filling. The protocol is inspired by ideas from dark state laser cooling, but 
replaces electronic states with motional levels, and spontaneous emission of photons by emission of phonons 
into a Bose-Einstein condensate, in which the lattice is immersed. In our model, achievable temperatures corre- 
spond to a small fraction of the Bloch band width, and are much lower than the reservoir temperature. 



PACS numbers: 03.75.Lm, 42.50.-p, 32.80.Pj 



5h 
43 



o 
o 



> 
in 

(N 

o 

O 

-i— > 



T3 

c 

o 
o 



X 



Fundamental advances in atomic physics are often linked to 
the development of novel cooling methods, as illustrated by 
laser and evaporative cooling, which led to the recent realiza- 
tion of degenerate Bose- and Fermi-gases 1 1 ]. This has further 
led to the achievement of strongly correlated atomic ensem- 
bles in the lowest Bloch band of an optical lattice 
However, in order to realise some of the most interesting con- 
densed matter phases predicted for lattice Hamiltonians, even 
better purification of the motional state is necessary, in partic- 
ular for atoms in a partially filled Bloch band |5]. Here we 
propose a method for cooling atoms to mean energies much 
smaller than the width of the lowest Bloch band 47°. In our 
setup (c.f. Fig. [U) lattice atoms a are excited to the first Bloch 
band via a Raman laser pulse except when they occupy Bloch 
states with quasi-momentum close to zero - so-called dark 
states. The lattice is immersed in a Bose-Einstein Condensate 
(BEC) of a different atomic species b, so that the atoms can 
subsequently decay back to the lowest band via collisional in- 
teractions with the BEC reservoir |6, 7]. Thus the atom is re- 
cycled to the lowest band by emission of a phonon - or more 
precisely, a Bogoliubov excitation [6]. By repeated applica- 
tion of laser excitation and "spontaneous emission", cooling 
into the dark state region of quasi-momenta is achieved with- 
out loss of atoms 11911 (c.f. Fig.^i). 

This method is inspired by the seminal Kasevich-Chu 
scheme |8] for sub-recoil laser cooling but replaces 

internal atomic states by Bloch band excitations, and sponta- 
neously emitted photons by phonons. Our scheme thus oper- 
ates on a much smaller energy scale than laser cooling, with 
correspondingly lower temperatures. This method can also 
be seen as a form of sympathetic cooling, where the energy 
is removed by phonons with energies equal to the Bloch band 
separation. Such phonon modes will initially be in the vacuum 
state, giving an effective T — reservoir, and allowing tem- 
peratures significantly lower than the BEC reservoir tempera- 
ture, in contrast to standard sympathetic cooling. The ability 
to switch the collisional interactions via Feshbach resonances 
lHUl enables us to study the cooling scenario in the weakly in- 
teracting gas, creating strongly correlated phases by ramping 
up the interaction in a final step. 

On a formal level our cooling scheme can be written as the 
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FIG. 1 : (a) Cooling setup: Atoms a in an optical lattice are coupled 
to the first excited motional state via a Raman process, and decay to 
the ground motional state due to collisional interactions with a BEC 
of species b in which the lattice is immersed. Tunnelling between 
neighbouring sites with amplitude J a gives rise to Bloch bands, (b) 
Momentum space picture: Atoms with higher quasi-momentum q 
are excited to the upper Bloch band, and decay to random quasi- 
momentum states. After several cycles, atoms a collect in a dark 
state region near q = with low excitation probability. 



iterative application of a map, 
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where the density operator p ; - describes the atoms in the low- 
est Bloch band before the jth step. Each step consists of 
two parts: the coherent laser excitation Ej of lattice atoms 
a, and the dissipative decay D returning atoms to the lowest 
band via coupling to the reservoir b (Fig.Q. To achieve the 
best possible cooling, differently shaped excitation pulses Ej, 
7 = 0,... ,N P — 1 are applied, and this sequence is repeated, 
with Ej = Ej mo( j w . This repeated application of the map cor- 
responds to a purification of the density operator, starting from 
an initial mixed state (e.g., a thermal distribution) towards a 
pure state (at zero temperature, T = 0). In order to find ap- 
propriate forms of the Raman pulses and the action of D, we 
analyse the dynamics of the lattice atoms and their interaction 
with the reservoir gas. 

We consider a one dimensional model for the motion of 
atoms a, which is readily generalised to higher dimensions. 
Including Raman coupling, the Hamiltonian is H a = Hq + Hj, 
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with 
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Here, A™ and (A™) T are annihilation and creation operators 
for quasi-momentum q in Bloch band a 6 {0, 1}, satisfying 
Bose or Fermi (anti-)commutation relations. The kinetic en- 
ergy is e™ = —2J a cos(qd), where d is the lattice spacing, J a 
are the tunnelling amplitudes (with J° > and J 1 < 0, see 
Fig- ID 5 )' an d 03 is the band separation. The effective Rabi fre- 
quency £2 = Q.r f dxexp(— i§qx)w l (x)w°(x), where £2# is the 
two photon Rabi frequency as a function of time during the 
pulse and w a (x) the Wannier functions. The Hamiltonian is a 
two-band model, written in a rotating frame with Raman de- 
tuning 8, and 5q denotes the momentum transfer. We choose 
units h = ks = l, where kg is the Boltzmann constant. The pa- 
rameters £2, 8 and Sq will change during the pulse sequence, 
but be constant during a given pulse j. Onsite interactions be- 
tween lattice atoms are represented by Hi, with hf the number 
operator for atoms in site / and band a, and U aa ' the associ- 
ated onsite energy shifts ||2|] . 

The density-density interaction between lattice atoms a and 
a three dimensional BEC reservoir b, which gives rise to the 
decay £>, is described by the Hamiltonian |6] 
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Here, the operator b\ creates a Bogoliubov excitation with 
momentum k = (k,k y ,k z ), and neglecting overlap of Wan- 
nier functions in different lattice sites, the coupling G„ , « 

gab(S(k)Pb/V) 1 / 2 Jd 3 xe iks w a (x)w a ' (x). The strength of the 
inter-species contact interaction is denoted g a b, pi, and V are 
the density and volume of the BEC reservoir, and 5(k) is the 
static structure factor 1 1]. For excitation energies less than the 
chemical potential /j, excitations are sound waves for which 
5(k) is strongly suppressed, and S(k) — > as |k| — > |1]. 
For energies larger than /j, excitations are in the particle-like 
sector of the spectrum, with much larger 5(k) — > 1. Here 
we will typically have 47° < /j < CO, so that decay between 
bands is induced by particle-like excitations with strong cou- 
pling, but collisional processes between the reservoir and 
atoms in the lowest Bloch band are suppressed. In close anal- 
ogy to Ref. 0, 0] we derive a master equation for the re- 
duced system density operator, describing the decay between 
bands in the Born-Markov approximation 1 20] . The associ- 
ated Liouvillian is L [p] = Y.k (?- c kPc\ — c^c^p — pc\ck\ /2. 

The momentum k along the lattice axis is bounded by \k\ < 
-\/2/M£,co due to energy conservation, where mt, is the mass 
of atoms b, and the jump operators Ck are defined as cj — 



^q^q-k(Ai^- ^ ne resu lti n g decay rates Tk for spontaneous 
emission of a phonon with momentum k projected on the 
axis of the lattice can be written explicitly for deep lat- 
tices, where CO 3> I-/ 1 ),/ and the individual lattice sites can 
be approximated as harmonic oscillator potentials. We find 
= gy h pt,m a aQk 2 exp(— a^k 2 /2) /2L, with ao the size of the 
ground state in each lattice site, m a the mass of atoms a, and 
L the length of the ID lattice. We denote the total decay rate 
from the excited band by r = J^I*. We consider a situation 
where dissipation is switched off during the excitation step, 
so that Ej and D occur separately. This can be achieved e.g., 
by tuning the collisional interaction so that g a i, ~ 0. We can 
read the action of Ej and D for a given step j from the master 
equation. 

We first illustrate the cooling process for a single lattice 
atom, designing a sequence of Raman laser pulses, where 
the j-th pulse excites the atom with initial quasi-momentum 
q in the lowest band to the first excited band with probabil- 
ity Pj{q). We require Pj(q) — for q w 0, but Pj(q) — > 1 
for states with high quasi-momentum (c.f. Fig. E 5 )- In 
analogy with Raman cooling schemes in free space |8] we 
choose square pulses with duration Xj — Tt/Q.,; for which 

Pj(q) = sin 2 (y^ +£l 2 x J /2)a 2 /(5 2 )+&ij +£22), with the 

effective detuning 8 q+Sq . = CO + E l q+Sq -E°-8j. 

An example of an efficient pulse sequence is shown in 
Fig. |2^-c. We begin with an intense laser pulse which res- 
onantly excites atoms with momentum q ~ n/d around the 
edges of the Brillouin zone (Fig. [2^). The subsequent pulses 
move the resonant transition closer to q = by adjusting the 
momentum transfer 8^ ; and Raman detuning 8 ; (Fig. |2j>). In 
order to prevent excitation of atoms with q — 0, we decrease 
£2 and increase the pulse duration, x for the later pulses, each 
time achieving Pj(q = 0) = 0. To resolve the transition we 
must always have £2 -C 8 {J 1 1 , and therefore x > 7t/8 1-/ 1 1 . Note 
that it is the value of J 1 and not J° that gives the resolu- 
tion of the excitation. However, the relationship between J 1 
and J° is fixed by the lattice depth (e.g., for the parameters 
used in Figs. |2jg| V = 10co fl , we have J° = 0.019co fi and 
/ = — 0.25<Mr). By combining a sequence of 5 pulses, one 
can efficiently excite most atoms with \q\ > 0, as shown in 
Fig-Et- Note that the pulse sequence has been carefully de- 
signed to avoid significant population outside the lowest two 
bands. 



To quantitatively analyze the cooling process we numeri- 
cally simulate the time evolution of the density operator p 
using quantum trajectories 11211 starting from a completely 
mixed state in the lowest band (T 47°), with ~ 10 5 trajecto- 
ries. During the cooling process, the momentum distribution 
develops a sharp peak near q = after very few iterations as 
shown in Fig. 01. In Fig. [2£ we plot the temperature of the 
state as defined by kgT /2 — ft 2 sin (Aqd)(J°) 2 /(OR, where Aq 
is the half- width of the momentum distribution at e~ 1 of the 
maximum value. We find excellent agreement between our 
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FIG. 2: (a)-(c) Excitation probability P — j(q) for a sequence of = 
5 pulses: first (a-c, dash-dot), second (b-c, dotted), and the remaining 
pulses (c). Parameters used: Q, = (27.9, 13.7, 13.7,2.37,2.37)7°, 
Sq = (0.16,-1.75, 1.75, -2.63, 2.63) /d, (8 - to) = 
-(28.4,25.8,25.8,24.7,24.7)7° for the different pulses, (d) Suc- 
cessive narrowing of the momentum distribution in the lower Bloch 
band after 0, 1 , 3, 5 and 10 cooling cycles from numerical simulations 
(M = 101 lattice sites) based on the pulse sequence in (a)-(c). We 
choose parameters for 87 Rb in the lattice and 23 Na in the reservoir, 
with r = 537° from m a jm\, = 3.73, p/, = 5 x 10 14 cm~ 3 and 
scattering length a a i, = 14nm. Vo = IOcor, and (Or = 2% x 3.8kHz. 
(e) Temperature vs. time for a single atom: crosses and circles 
denote numerical, solid lines analytical results based on Levy 
statistics. Pulse sequences for circles: same as in (d); crosses: 
N p = 3 pulses with £^ = (32.6,7.9,7.9), Sqd = (0.31,2.12,-2.12), 
(8 -co) =-(28.4,25.3,25.3)7°. 



numerical calculations and analytical results obtained with 
Levy statistics as a function of time. The latter pre- 

dicts a final temperature T <*t~ l for square pulses, as shown in 
Fig-Et- For a zero-temperature reservoir, and parameters as in 
the caption for Fig.|2] we reach temperatures 7~2x 10 _3 7° 
in time t f J° — 30. 

Finite temperature Tf, in the reservoir can lead to sympa- 
thetic heating of lattice atoms a by absorption of thermal 
phonons, as described by Hi nt . However, this process is 
forbidden by energy and momentum conservation, provided 
J° < \J fi(i3Rm a / (2nib) /Jt. In detail, energy conservation re- 
quires c|k| = e° — £°, and conservation of momentum along 
the lattice direction leads to k = q — q' {\k\ < |k|), where the 
atom a is scattered from quasi-momentum q as — > q' by ab- 
sorption of a phonon with momentum k, and c is the sound ve- 
locity in the BEC. These conditions cannot be fulfilled unless 
the above inequality is violated. Higher order processes in- 
volving two or more thermal phonons will be small. In numer- 
ical simulations we also checked that the cooling protocol is 
insensitive to small timing errors. Whilst in the above protocol 
we have switched off the decay during application of Ej, we 
can leave decay switched on, provided that T <C 1/x <SC \ J l \. 
This will restrict the length of the pulses that can be applied, 
thus slowing the cooling process. 



The cooling scheme can be readily adapted to many bosons 
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FIG. 3: Numerical simulation of the QBME: (a) Temperature as 
a function of cooling time for bosons (crosses) and fermions (cir- 
cles), (b) Momentum distribution in band a = for fermions after 
(dashed), 1 (dash-dot), 2 (dotted), and 20 (solid) cooling cycles, 
each with N p = A pulses. Parameters used: Bosons: As for Fig|2^-c, 
but A' = 51 particles. Fermions: N = 71, M = 101, Blackman pulses 
withxT = (1.78,1.78,6.8,6.8), Sqd = (0.19,-0.19,0.75,-0.75), 
(co- 8) = (28.4,28.4,27.9,27.9)7°; V = 10co R , co R = 2% x 6kHz, 
m a /m b = 1.74 and T = 52.67°. 



or fermions. For bosons, we assume that the collisional in- 
teraction between atoms a is tuned to zero (Hj — > 0). We 
work out the efficiency of the cooling protocol by deriv- 
ing a quantum Boltzmann master equation (QBME) fl3ll . 
which describes transitions between classical configurations 
of atoms occupying momentum states in the Bloch bands, 
m = [{m^}g,{m^} ? ], where m" is the occupation number of 
quasi-momentum state q in band a. This is derived from the 
master equation by projection of the density operator p onto 
diagonal elements, £PpiP = £ m w m |m) (m|, neglecting off- 
diagonal coherences. For the excitation step, Ej, the evolution 
is computed exactly from the excitation probability Pj(q), and 
for the decay, D, we obtain 



k.q 



where m' = m — v q -k,q is the resulting configuration when a 
particle with quasi-momentum q in the upper band decays to 
the lower band with new quasi-momentum q — k, i.e., e q -k.q 
is a configuration vector with tn® k = 1, m f ' = — 1 and all 
other entries zero. The upper (lower) signs are for bosons 
(fermions). The approximation inherent in neglecting off- 
diagonal coherences only plays a role during the decay step, 
where these coherences couple to the occupation probabilities. 
We remark that an exact physical realisation of the QBME can 
be obtained by modulating the lattice depth after each excita- 
tion step, randomising the off-diagonal elements Jl4ll . 

Fig.[3k shows the decrease in temperature as a function of 
time for bosons and fermions, obtained from monte carlo sim- 
ulations of the QBME II 311 . For bosons, we use the same exci- 
tation pulse sequence as for a single atom in Fig.|2] The cool- 
ing process in this case outperforms that for a single atom, 
reaching low temperatures on shorter times due to bosonic 
enhancement (here we compute temperature as for a single 
atom, but fitting a Gaussian to the q ^ momentum distri- 
bution). For fermions, the pulse sequence must be changed 
to create a dark state region of quasi-momenta with \q\ < q?, 
where qp is the Fermi momentum, in order to cool towards a 
T = Fermi distribution. In this case, time-square pulses are 
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FIG. 4: (a) Momentum distribution in the lowest Bloch band from 
numerical simulations after a single excitation pulse beginning with 
equal occupation, for interaction strengths U aa = (0,0.1,0.3)7° 
(solid, dashed, and dotted lines). Parameters used: N = 5 atoms 
in M = 41 lattice sites 8qd = 1, £2 = 1.057°, (8 - co) = 27.97°, 
V0 = 10c0s, COfl = 271 x 3.8kHz. (b) Energy £ t| of the state obtained 
by beginning in the U = ground state with N = 10, M = 21, and 
ramping the interaction strength as from U 00 ~ to U ~ 207° in 



time l r , as U 



oo . 



:207°(l-{l+exp[(f-T r /2)/(t r /10)]}- 1 ). 



no longer efficient as there is a large secondary peak in Pj(q) 
(see Fig.lJJ}), and we instead use Blackman pulses |8], which 
approach Pj(q) = monotonically. The momentum distribu- 
tion develops the expected shape of a cold Fermi distribution 
after very few iterations (Fig.|5J)) (we compute temperatures 
for fermions by fitting a Fermi distribution to these results). 

Our model predicts that the temperature will always de- 
crease with increasing cooling time. Experimental imperfec- 
tions will, in practice, give rise to decoherence and heating 
(e.g., from spontaneous emissions |2], collisions with back- 
ground gases, scattering multiple phonons). One assumption 
made above for bosons was that the interaction between lat- 
tice atoms a is approximately zero. Using time dependent 
DMRG methods 1 15] we computed the population remaining 
in the lowest band after an excitation pulse with interactions 
present (Fig.|4k). For U aa <C 1/t <C I7 1 ! there is no signifi- 
cant change in the excitation profile and the above conclusions 
should not change. 

Finally, after the cooling we adiabatically ramp up the in- 
teraction strength to produce an interacting system, assuming 
that the system is decoupled from the reservoir b. This is illus- 
trated in Fig.|4]for the case of ramping from a non-interacting 
to a hard-core Bose lattice gas in ID (Tonks gas). We compute 
the evolution from a Bose-Hubbard model (H a with no Raman 
coupling, and only the lowest band), again using time depen- 
dent DMRG methods 1131 . The energy deposited during the 
ramp of the interaction strength is plotted as a function of the 
ramp time x r , and for ramping times on the order of 10/7°, we 
observe negligible heating within the Bloch band. Numerical 
tests using example excited states also showed that the energy 
difference between different initial states is not significantly 
increased during the ramp. 



In summary, filtering quasi-momentum states in the low- 
est Bloch band of an optical lattice and recycling them by 
"spontaneous emission" of phonons combine to give a cooling 
scheme producing temperatures a small fraction of the lowest 
Bloch band width. These temperatures are a necessary step 
towards the realisation of strongly correlated quantum phases 
not currently accessible in optical lattice experiments. 
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